Correlation effects on a topological insulator at finite temperatures 
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We analyze the effects of the local Coulomb interaction on a topological band insulator (TBI) by 
applying the dynamical mean field theory to a generalized Bernevig-Hughes-Zhang model having 
electron correlations. It is elucidated how the correlation effects modify electronic properties in the 
TBI phase at finite temperatures. In particular, the band inversion character of the TBI inevitably 
leads to the large reduction of the spectral gap via the renormalization effect, which results in the 
strong temperature-dependence of the spin Hall conductivity. We clarify that a quantum phase 
transition from the TBI to a trivial Mott insulator, if it is nonmagnetic, is of first order with 
a hysteresis. This is confirmed via the interaction dependence of the double occupancy and the 
spectral function. A magnetic instability is also addressed. All these results imply that the spectral 
gap does not close at the transition. 
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I. INTRODUCTION 

Recently, the topological band insulator (TBI) has re- 
ceived much attention as a new research platform in 
condensed matter physics [l|, The TBI has a non- 
trivial band structure resulting from spin-orbit interac- 
tions (SOIs). A remarkable property of this phase is that 
while it has the charge excitation gap in the bulk, an odd 
number of gaplcss edge states appear at boundaries Q, 
which are robust against nonmagnetic impurities due to 
time reversal symmetry. The non-trivial band structure 
in the TBI was originally proposed foi_graphene and 
also for HgTe/CdTe quantum wells Q. The latter TBI 
was confirmed by transport measurements [lol |. After- 
wards, TBIs in three dimensions were proposed theoreti- 
cally [7|-|9| , and confirmed in several bismuth-based com- 
pounds through the angle-resolved photo-emission spec- 
troscopy measurements which elucidated the existence of 
the gapless edge states [lTI - [l3j . 

While the TBI can be described by a non-interacting 
band theory, correlation effects on this phase have at- 
tracted much attention in these days, since novel aspects 
of electron correlations would appear under topologically 
nontrivial conditions. This issue has further been stimu- 
lated by the fact that there are some iridium oxides, such 
as Na2lr03 [11] and AalraOy (A = Pr,Eu) [Hi, which 
could be candidates of TBIs with strong electron cor- 
relation. Theoretical studies on the correlation effects 
done thus far are roughly categorized into two types. 
One is concerned with the case where the nonlocal re- 
pulsive interaction could induce a quantum phase tran- 
sition, thus leading to topological insulators, as demon- 
strated theoretically for both 2D and 3D systems with 
mean field analysis [l6l - [l9| . The other focuses on how 
the TBI phase is changed under strong electron correla- 
tions. For example, a competition between an antiferro- 
magnetic (AF) phase and the TBI phase was studied by 
numerical approaches [20| - (23 | as well as mean field the- 
ory [l^l- According to a quantum Monte Carlo study by 
Hohcnadler et al. [2l| , it was clarified that the TBI phase 



can change into a topologically trivial AF phase, which 
was supported by the bulk phase diagram obtained by a 
variational cluster approach and also a cluster dynamical 
mean field study [2j, |2J] . The former study also clarified 
the absence of edge states in the AF phase [2^ . A com- 
petition of the topological phase and a charge-density- 
wave phase was also studied [2^, [13] . According to the 
exact diagonalization analysis, the nearest neighbor in- 
teraction induces a first-order transition, and the TBI 
phase changes into the charge-density-wave phase. 

In spite of these intensive studies, it is still not clear 
how the electron correlations affect the behavior of the 
physical quantities characterizing TBIs such as the spin 
Hall conductivity, since most of the previous studies have 
focused on the bulk gap and the edge states at zero tem- 
perature. This issue is particularly important to discuss 
finite-temperature properties because the strong correla- 
tion should largely change the relevant energy scales of 
the system, bringing about a strong renormalization of 
the physical quantities. 

The purpose of this paper is twofold. We first elucidate 
how the correlation effects modify electronic properties in 
the TBI phase at finite temperatures with particular em- 
phasis on the spin Hall conductivity. The deviation of the 
spin Hall conductivity from the quantized value at finite 
temperatures is clarified in terms of electron correlations. 
To this end, we employ a generalized Bcrncvig-Hughes- 
Zhang (BHZ) model having local electron correlations, 
and treat it with dynamical mean field theory (DMFT) 
[28l - [30j . It will be revealed that the spectral gap of the 
TBI is gencrically reduced by the correlation effect, as 
a result of band inversion, leading to the strong tem- 
perature dependence of the spin Hall conductivity. We 
then discuss how the TBI changes into the MI that could 
emerge in a strongly correlated region. We will show 
that a first-order transition occurs between the TBI and 
the Mott insulator (MI) , as far as the nonmagnetic Mott 
phase is concerned, and thus a gap closing does not oc- 
cur at the TBI-Ml transition. The magnetic instability 
is also discussed. 
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II. MODEL AND METHOD 



from the following equations |31|-|3 



We extend the BHZ model to include on-site Coulomb 
interactions. The Hamiltonian reads 
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where n^^a^o. = cj The operator cl^ ,^{ci^a,a) 
creates (annihilates) an electron at site i and orbital 
a = 1,2 in spin state a =t, 4-- Off-diagonal elements 
of the hopping matrix t drive the system into a non- 
trivial band insulator in the non-interacting case. The 
angle 9 in the hopping matrix t specifies the direction 
of hopping (see Fig. [Ha)). We analyze the system 
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where Z/oc is the partition function for the effective 
impurity model. In this method, Eq. ([6]) is expanded 
with respect to the hybridization functions i^a.a' j and 
the resulting integral is evaluated with the Monte 
Carlo method. Absence of the Trotter decomposition 
makes it easier to access the low temperature region. 
For simplicity, we study the particle-hole symmetric 
case and choose the model parameters as ii = t2 = t, 
tsp = 0.5i and £1(62) = The hopping integral t is 



chosen as the energy unit. 



III. RESULTS 
A. spin Hall conductivity 




FIG. 1: (a) Sketch of the BHZ model on a square lattice. The 
hopping integral between orbitals depends on 9 which defines 
the angle between hopping direction and x-axis. (b) Sketch 
of the DMFT. The lattice model is mapped to an eflfective 
impurity model. 

with DMFT, which treats local correlations exactly and 
is suitable for the systematic calculation for arbitrary 
strength of the Coulomb interaction. In DMFT, the orig- 
inal lattice problem is mapped onto an effective impu- 
rity model, which is solved self-consistcntly. The self- 
consistent equation for a paramagnetic phase is given by 

where /i(j(k) is the Fourier transform of the hopping 
matrix. The self-energy of the lattice Green's function 
E(j(ia;) can be computed from the Green's function g{iLu) 
for the effective impurity model. 

To solve the impurity problem, the continuous-time 
quantum Monte Carlo method (CT-QMC) is employed 
here. In this method, the full Green's function is obtained 
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FIG. 2: (Color online) The spin Hall conductivity as a func- 
tion of the interaction strength. Solid circles, open triangles, 
and open circles denote the conductivity at T = 0.04, 0.05, 
and 0.08 respectively. Inset: spin Hall conductivity as a func- 
tion of temperature. Open circles show the data at [/ = 8. 
The solid red line shows the data calculated with the renor- 
malized parameters which are obtained by low-energy expan- 
sion at (7 = 8. 

Let us first discuss the spin Hall conductivity, which 
characterizes the topological nature of the system. The 
spin Hall conductivity is calculated with the following 
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formula, 



xy 
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where V denotes the volume of the system. Here, some 
comments are in order on the calculation of the conduc- 
tivity. As well as in the non-interacting case, the con- 
ductivity should be a multiple of /h at zero temper- 
ature, which is checked as follows. Taking into account 
that the conductivity is antisymmetric tensor and that 
the vertex part is connected through the Ward identity 
(A^(p,p) = g^G~"'^(p)), we may represent the conduc- 
tivity as the following equation: 
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where the notation p = (ia;,p) is used. N is the first 
Chern number expressed in terms of the Green's func- 
tion [131 ■ Therefore, even in the correlated system, the 
conductivity should be a multiple of e'^/h at zero tem- 
perature [H, [H]. Besides, in the DMFT framework, the 
self-energy has no momentum dependence, and we can 
replace the vertex part with dh^/dk^. Thus, the con- 
ductivity, calculated via eq. (|10l) . is exact within the 
DMFT framework. The analytic continuation is approx- 
imately done by replacing Im dK^{uj + iS)/duj with Im 
dK{iLu)/diuj. 

Note that in the above expression of erf , there are no 
current vertex corrections from quasiparticlc interactions 
which are related to the backflow. and normally exist in 
the expressions of transport coefficients. This is because 
that the current vertex corrections are associated with 
current dissipations raised by quasiparticlc scatterings, 
and thus they do not appear in the expression of the 
dissipationless Hall conductivity with which we are con- 
cerned here. Therefore, the above formula for the spin 
Hall conductivity is exact even for interacting electron 
systems. 

We now turn to the numerical results. The calculated 
spin Hall conductivity is plotted in Fig. [5] as a func- 
tion of the interaction strength at several temperatures. 
It is seen that the conductivity in the non-interacting 
case takes values close to unity in units of /h with lit- 
tle temperature dependence, implying that the system is 
indeed in the TBI phase. When the interaction is intro- 
duced, the conductivity gradually decreases and shows 



a strong temperature dependence with increasing U. It 
eventually vanishes in the large interaction region. Note 
that the conductivity in the TBI region is expected to 
reach unity at zero temperature. Actually, as seen in 
the inset of this figure, the conductivity for U/t = 8 
(open triangles) increases toward unity with decreasing 
temperature. We can confirm this tendency by focus- 
ing on the low-energy properties of the self-energy due 
to the Coulomb interaction; the computed values exhibit 
the almost identical behavior to the solid red line calcu- 
lated with the renormalized parameters, which are ob- 
tained by the low-energy expansion of the self-energy as 

E^(a;) = S^(0) + ^(O)a; + 0{uj^). At T < 0.05, no 
strong temperature dependence of the renormalized pa- 
rameters is observed. We thus conclude that the spin Hall 
conductivity should reach unity at zero temperature and 
therefore the TBI should be stable for U/t < 12. 

It is instructive to point out that the above character- 
istic temperature/interaction dependence of the conduc- 
tivity is not specific to the model parameters employed 
here, but more generic for correlated TBIs. As will be 
discussed momentarily below, the renormalization of the 
charge gap occurs due to interactions, which may increase 
the effective temperature with respect to the renormal- 
ized gap. This tendency is expected for the TBIs, so that 
the decrease of the conductivity with increasing temper- 
ature/interaction should be generally observed in TBIs. 

In order to clarify the origin of the strong-temperature 
dependence of the conductivity, let us now look at the 
one-particle spectral function. The spectral function 
(A(w) = -^lmTTG{u} + i5)) is calculated at T/t ^ 0.05. 
The analytic continuation procedure to convert the Mat- 
subara frequency to the real frequency is done numer- 
ically with the maximum entropy method (MEM). The 
results obtained for several choices of interaction strength 
U are plotted in the upper panel of Fig. |3l In this figure, 
we can see the reduction of the gap and the evolution of 
Hubbard peaks with increasing [/, which is consistent 
with the previous studies on the Kane-Mele-Hubbard 
model [111, m, mi • In the large U region, the renormal- 
ized peaks in the low-energy region eventually disappear, 
implying that a quantum phase transition occurs between 
the TBI and the MI. 

We note here that in generic band insulators the spec- 
tral gap is not necessarily reduced via the renormaliza- 
tion. Sometimes it is even enhanced by the interaction. 
Actually, according to Ref. [s^l , if the system has a local 
hybridization, the off-diagonal self-energy increases the 
local hybridization and hence enhances the gap. In this 
sense, whether the gap is enhanced or reduced depends 
on the model systems employed. However, a remarkable 
point is that the gap for the TBIs should be always renor- 
malized to be smaller due to electron correlations. This 
intriguing property is a consequence of the band inversion 
inherent in the TBIs. Note that the TBIs have the gap 
resulting from the SOIs. If the interaction is switched 
on, the SOIs should be renormalized because the SOIs 
contribute to the kinetic energy. This renormalization 
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FIG. 3: (Color online) Upper panel: Local density of states for 
each value of U. Figs, (a), (b) and (c) correspond to the TBI 
phase, while Figs, (d), (e) and (f) correspond to the MI phase, 
respectively. Note that in 11.5 <U/t < 12.125, a hysteresis is 
found. Lower panel: Renormalized energy sphtting between 
the two orbitals. 
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FIG. 4: (Color online) Double occupancy of each orbital as 
a function of interaction: (a) T/t^0.05 (b) T/t=0.12. Note 
that at r/t=0.12, no hysteresis is found. 



gives the reduction of the spectral gap, as shown in Fig. 
[21 Although the renormalization of the energy sphtting 
62 — e*, shown in the lower panel of FiglJl seems to cause 
the gap renormalization, this renormalization should be 
irrelevant since it is essentially absorbed in the energy 
shift due to the band inversion in the TBIs. Therefore, 
we can say that the reduction of the gap due to the inter- 
action and the resulting strong-temperature dependence 



of the spin Hall conductivity are interesting generic fea- 
tures common to the TBIs. This result also implies that 
strong electron correlation is harmful to the topological 
stability of the quantum spin Hall effect. 



B. Mott transition 

We now discuss what happens in the strong correlation 
region. In Fig. [31 we have already observed the disap- 
pearance of the renormalized peak structure for the TBI 
in the spectral function, but the nature of the transition 
is not clear yet. To address this issue, we compute the 
double occupancy of electrons at each site, which is a key 
quantity to characterize the Mott transition. In Fig. [H 
the double occupancy of each orbital is plotted as a func- 
tion of the interaction strength. It is clearly seen in Fig. 
[31 (a) that the double occupancy exhibits a jump and a 
hysteresis around U/t ^ 12. This behavior indicates that 
the TBI changes into the MI via a first-order transition. 
Therefore, in the upper panel of Fig. [31 the coherence 
peaks in the low energy region vanish suddenly, and the 
system acquires the Hubbard gap in the MI phase. This 
is contrasted to the non-interacting case, where the gap 
closing occurs continuously when the TBI changes into 
a trivial insulator. In the MI phase, the gap size is of 
the order of J7, so that the temperature effect is not im- 
portant in this region. As seen in Fig. [2l the spin Hall 
conductivity is almost zero, which confirms that the MI 
belongs to a trivial phase. 

In Fig. m (b), the results at higher temperatures are 
shown. It is seen that the first-order transition is now 
changed into a crossover where the jump and hysteresis 
are smeared. Note that this kind of crossover behavior is 
usually accompanied by the first-order Mott transitions 
at higher temperatures. 
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FIG. 5: (a): Sketch of the spin configuration in the AF phase. 
Electron spins are denoted by the blue arrows. The numbers 
1, 2 specify the orbital indices, (b): Magnetic moment as a 
function of the interaction at T/t = 0.05. 



C. magnetic instability 

We have so far restricted ourselves to a paramagnetic 
MI phase. Note that in the DMFT framework, the as- 



sumption of the paramagnetic MI phase corresponds i 
imphcit introduction of frustration, which is not real 
included in the present model (see the discussions b 
low). Now, let us turn to the magnetic instability of tl 
system. From the estimation of the exchange intera 
tions with a perturbation theory in hopping, it is four 
that the off-diagonal elements of the hopping matrix 
provide the exchange interactions which depend on sp 
components; the z-component of the exchange intera 
tion is antiferromagnetic (AF), while the in- plane cor 
ponents are ferromagnetic (AF) between neighbors in ; 
(x-) direction. This results from the phase factor 9 ar 
spin dependency of the hopping matrix. From the abo^ 
considerations, the spin configuration of the AF pha; 
is expected to be the one shown in Fig. [S] This Ion 
range AF ordered phase is treated with the sublattii 
method in the framework of DMFT. The order param 
ter, which is defined by to = X]a("'ct.t ~ 'T'Q,i)/2, is plo 
ted in FiglSj^b). In this figure, we can find a second-ord 
quantum phase transition. By performing similar calc 
lations for various sets of parameters, we end up wi1 
the phase diagram shown in Fig. [HI It is seen that tl 
paramagnetic Mott transition is now masked by the A 
phase. This kind of behavior is commonly observed : 
the DMFT treatment for the MI transition. Therefor 
we naturally expect the transition from the TBI to tl 
AF insulator for our generalized BHZ model. The abo^ 
results imply that the spectral gap which persists up i 
[/ ~ 12t in the paramagnetic case docs not close at tl 
magnetic transition aX U 5t. This behavior is simili 
to the result obtained for the Kane-Mele- Hubbard mod 
in Ref.jllj, though another research work reports thi 
the gap closing occurs at the transition [2^. 

If the system is geometrically frustrated, however, tl 
paramagnetic MI phase mentioned in the previous su' 
section might be stabilized. Here, we make a comment c 
whether we can have a chance to expect the Mott tra: 
sition in such frustrated systems. A frustrated system : 
the BHZ model can be realized if the next nearest neig' 
bor (NNN) hopping is taken into account in the square 
lattice. We have confirmed that in the particle-hole sym- 
metric case, the system is in the TBI phase a.t U = if 
the NNN hopping meets the condition: £2 — 61 < A{t — t'), 
where t' represents the NNN hopping. Therefore, we 
could have a quantum phase transition between the TBI 
and the nonmagnetic MI discussed here, starting from 
the frustrated BHZ model. This problem is now under 
consideration. 

In this connection, we wish to mention that there are 
some other candidates for TBIs which have geometri- 
cally frustrated lattice structure. A certain extension 
of a kagome lattice, which has a non-trivial topological 
band structure [39j, and pyr ochlore iridium oxides such 
as A2lr207 (A = Pr,Eu) \m. [ssj are proposed as candi- 
dates of frustrated TBIs. In these frustrated systems, we 
may expect that the first-order quantum phase transition 
should occur between the TBI and the MI, if we could 
realize the TBI in the weak coupling region. 
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FIG. 6: (Color online) The temperature vs. interaction phase 
diagram. The dotted blue line denotes the second order tran- 
sition between the TBI phase and the AF phase. If the AF 
phase is suppressed, a first-order Mott transition could be 
observed (see text). Solid red lines denote the Mott transi- 
tion line, where the coexistence phase accompanied by the 
hysteresis is found in the region surrounded by red lines. 



Before closing this section, we make a brief comment 
on a possible topological aspect of the antiferromagnetic 
phase. It was previously proposed that a topologically 
nontrivial phase could appear in the antiferromagnetic 
region [40| . However, we have not found the non-zero 
spin Hall conductivity in this phase within the present 
approach. Detailed study on this issue is interesting, 
which is to be done in the future work. 



IV. SUMMARY 

In this paper, we have studied the effects of the lo- 
cal Coulomb interaction on the TBI by applying DMFT 
and CT-QMC to a generahzed BHZ model. We have 
calculated the spin Hall conductivity and found that the 
strong temperature dependence of the spin-Hall conduc- 
tivity results from the renormalization of the gap due to 
electron correlations. We have pointed out that the re- 
duction of the gap via the renormalization is a generic 
feature common to the TBIs, and is raised by a band 
inversion structure characteristic of TBIs. We have also 
shown that the TBI-MI transition is of first-order if we 
restrict ourselves to a paramagnetic Mott phase. This 
condition corresponds to taking into account frustration 
in the DMFT framework. If the magnetic instability is 
considered in our generalized BHZ model, the Mott tran- 
sition is masked and the transition becomes of second or- 
der. Our results imply that a gap closing process does 
not occur for both cases of the paramagnetic Mott tran- 
sition and and the magnetic transition. We have also 
discussed possible situations in which the paramagnetic 
Mott transition would be observed for the systems with 
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frustration. It is an interesting open issue to treat the 
Mott transition for such frustrated systems in topologi- 
caUy nontrivial conditions, such as Iridium oxides. The 
appUcation of DMFT calculation to such systems is now 
under consideration. 
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